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Recent experiments on atom loss in ultra-cold Fermi gases all show a maximum at a magnetic 
field below Feshbach resonance, where the s-wave scattering length is large (close to inter-particle 
distance) and positive. These experiments have been performed over a wide range of conditions, 
with temperatures and trap depths spanning over three decades. Different groups have come up with 
different explanations, among them the emergence of Stoner ferromagnetism. Here, we show that 
this maximum is a consequence of two major steps. The first is the establishment of a population of 
shallow dimers, which is the combined effect of dimer formation through three-body recombination, 
and the dissociation of shallow dimers back to atoms through collisions. The dissociation process 
will be temperature dependent, and is affected by Pauli blocking at low temperatures. The second 
is the relaxation of shallow dimers into tightly bound dimers through atom-dimer and dimer-dimer 
collisions. We have constructed a simple set of rate equations describing these processes. Remark- 
ably, even with only a few parameters, these equations reproduce the loss rate observed in all recent 
experiments, despite their widely different experimental conditions. Our studies show that the lo- 
cation of the maximum loss rate depends crucially on experimental parameters such as trap depth 
and temperature. These extrinsic characters show that this maximum is not a reliable probe of the 
nature of the underlying quantum states. The physics of our equations also explains some general 
trends found in current experiments. 



I. INTRODUCTION 



Three-body processes occupy a unique place in the study of ultra-cold atomic gases. On the one hand, while 
two-body elastic collisions are important for thermal equilibrium, three-body processes are crucial for producing 
molecules of interested, such as those in the BEC-BCS crossover. On the other hand, three-body processes can also 
lead to atom loss from the trap, which limits the lifetime of the sample and constraints the time scale over which 
experiments can be performed. In addition, three-body process can also be used as diagnostic tool. Recent advances 
in Efimov physics are resulted from using three-body loss to identify the presence of trimers in a bosonic gas of 7 Li[TJ[5]. 

Theoretically, the study of three-body problem dates back to the early days of quantum mechanics and has 
revealed many exciting new phenomena, such as Thomas effect [3] and Efimov physics [3]. In the case of interest 
here, namely, an equal population of two hyperfine states of fermionic quantum gas such as 6 Li, Efimov physics is 
irrelevant. What is important is the process of three-body recombination, in which three atoms collide to form a 
shallow dimer, with the third atom carrying away the energy released in the process. We shall refer this as Process 
(I) and denote it as + + — >• D + A+ , and A^ + A^ + A^ — >• D + Ai, where D is the shallow dimer; A^ and 
A^ are the atoms in different hyperfine states which we simply denote as f and J,. Note that process (I) only occurs 
when the two-body s-wave scattering length a s is positive, since bound states only exist for a s > 0. The reverse of 
Process (I) is the dissociation of shallow dimers back to atoms, A^ + D — >• A^ + A^ + A±, which we refer to as Process 
(I'). The combined effect of (I) and (I') is to provide a net population of shallow dimers. 

Once shallow dimers D are formed, they can further relax into deep bound states in the van der Waals potential 
through either atom-dimer collision, A + D — > D* + A (referred to as process (II)), or dimer-dimer collision, 
D + D — >• D* + D* (referred to as process (III)), where D* is the deep bound state. These two processes are the key 
reasons for particle loss in a trapped gas, since both types of particles (dimer and atoms) in the final product will 
have a very high kinetic energies due to energy-momentum conservation, and will leave the trap. It should be noted 
that these two processes can only be activated after a population of shall dimers is formed as a result of Process 
(I) and (I'). Our study shows that the atom loss rate depend on a variety of factors such as trap depth and the 
temperature of the system. The location of maximum loss rate is not "intrinsic" in the sense that it can change by 
varying external conditions such as trap depth, even though the gas parameter fcpa s is kept fixed, where kp is the 
Fermi wave vector in the center of the trap. 



The paper is organized as the following. In Sec 



TTJ we give an overview of the recent experiments and the range 



of parameters appeared in these studies. In SecjHIl we discuss the processes (I) to (III) mentioned above. We show 
that on general grounds, the loss rate will have a maximum at a magnetic field below Feshbach resonance. In Sec |IV| 
we implement these processes into a set of rate equations, and give explicit expressions for the rate constants in these 
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equations. In Sec|Vj we show that the results obtained from the rate equations provide a good description for all 
current experiments which are performed under a variety of external conditions. In the concluding Sec VI we further 
discuss the implication of the "two step" process (i.e. [(I) + (I')] and [(II) + (III)] ) mentioned above. 



II. OVERVIEW OF RECENT EXPERIMENTS 



The experiments we are going to discuss are performed in a mixture of two-component Fermi gases with density 
n interacting with an s-wave scattering length a s . Interaction strength of the system is characterized by the 
"gas parameter" kpa s , where kp = (S^n) 1 / 3 is the Fermi wavevector. For positive scattering length, a s > 0, a 
two-body system can be either in a scattering state or in a bound state. In the literature, these bound states are 
sometimes referred to as "Feshbach molecules", or sometimes simply "dimers". The binding energy of the dimer 
is Ef, = h 2 /ma 2 s > 0. When the two-body system is in the scattering state (or in the bound state), it is referred 
to as in the "upper branch" (or "lower branch"). The generalization of this terminology to the many-body case 
(for a s > 0) is the following. If the quantum gas is in thermal equilibrium, it is referred to as in the "lower" 
branch. If the quantum gas is prepared in a state where dimers are absent, such as the case at sufficiently high 
temperature or sufficiently small kpa Sl it is referred to as in the "upper branch". By definition, the system in the 
upper branch is not in its ground state. However, if the rate of production of Feshbach molecules is sufficiently low, 
the atoms in the system can be regarded in a quasi-equilibrium state during the time when experiments are performed. 

To our knowledge, all recent experiments on atom loss in a trap are performed in the upper branch. These 
experiments include many studies of Grimm's group at Innsbruck |10j : the work of Jin's group at JILA [7j, the work 
of Ketterle's group at MIT [HE], and Salomon's group at ENS[5). Those experiments were performed under a wide 
range of conditions, with temperature T spanning over two orders of magnitude, trapping depth Vq over three orders 
of magnitude, and a wide range of Fermi temperatures Tp. The parameters of these experiments are collected in Table 
HJ In all these experiments, a maximum loss rate is found at a magnetic field B* below the magnetic field B^ where 
Feshbach resonance occurs. The positive scattering length at which the maximum loss takes place will be denoted as 
a* . The binding energy of the Feshbach molecule (or dimer) at a* will be denoted as E£ . 



TABLE I: Experimental Parameters. 



VojfiK) T(fiK) T F {/iK) k F a* s B* (G) ggj EZ(jj.K) E£/T jg/gj Vb/T Ref 



7.1 


0.28 


1.4 


2 


790 


4318 


0.7 


2.5 


0.098 


25.35 


MIT [5], Li 


10 


0.67 


0.79 


0.83 


201.4 


990 


1.64 


2.44 


0.164 


14.9 


JILA[7], K 


175 


22 


21 


1.6 


680 


700 


17.88 


0.81 


0.102 


7.95 


MIT [6], Li 


350 


22 


2.8 


0.359 


644 


427 


64 


2.9 


0.183 


15.9 


Innsbruck|10j. Li 


500 


30 


2.8 


0.319 


636 


380 


83 


2.76 


0.166 


16.6 


Innsbruck|10). Li 


1000 


60 


2.8 


0.286 


629 


342 


104 


1.73 


0.104 


16.6 


Innsbruck|10j. Li 


a 


2.4 


6 


1.2 


720 


1207 


7 


2.9 






ENS [8], Li 



"Ref. [8] stated that the trap depth is of the same order as temperature. The specific value for it was not reported. 



Different groups have different views on the physical origin of the maximal loss. The Innsbruck group interpreted 
their data as due to a two-step process, similar to the mechanism presented below, but without detailed formulation 
and calculations [TU]. We shall discuss this picture in detail below. The JILA group noted that at the maximal 
loss, the heating of the system also reaches a maximum. However, no specific proposals were made for the origin 
of the observation. The ENS group considered the maximum loss as a consequence of the rising three-body 
recombination rate at small a s , and the decreasing binding energy of the Feshbach molecule as a s approaches infinity 
at resonance. It is argued that the latter reduces heating (and hence loss rate) because the energy release in the 
3-body recombination is small. However, no quantitative comparison had been made with experimental data using 
this picture. In the earlier experiment of the MIT group, no explanation was given to this maximum loss. In 
ref. [5], which is performed at very low temperatures, the maximum loss rate is considered as evidence for Stoner 
ferromagnetism. Their point is that if the system turns ferromagnetic, different spins will segregate, which will 
naturally leads to a vanishing three body loss rate. The only problem, however, is that no segregation of up and 
down spins have ever been observed. This explanation is very different from all previous pictures, for it makes use of 
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intrinsic property of the ground state rather than specific microscopic scattering processes. It raises the general ques- 
tion of how reliable it is to use the 3-body loss as a tool to probe of the nature of the ground state, ferromagnetic or not. 



III. THE TWO-STEP PROCESS 

Let us first examine the loss channels available for the trapped gas. For simplicity, we start with a system consists of 
only atoms, namely, no shallow dimers are present. In that case, the only loss channel at time t = is the three-body 
recombination (process (I)), which we schematically denoted as 

Af + Af + Ai ->• Af + D Process (I) (1) 

and a similar process where up and down spins are interchanged. Here, A^(A^) stands for atom with spin up 
(down) atoms and D is a shallow dimer consists of two opposite spins. We have only considered those three-body 
recombinations leading to shallow dimers. In principle, three-body recombinations can also lead to deep bound states 
that exist in the two-body van der Waals potential. However, the Franck-Condon factors for these transitions are 
so much smaller than those for the shallow dimers that they render the transitions to deep bound state negligible 
in comparison. Process (I) is an exothermic process. An important point to note is that the energy released in the 
process is of the order the dimer binding energy Ef, which is typically much smaller than the trap depth Vq, i.e. 
Eb <C Vb (see Table [if. As a result, the shallow dimers (D) formed in process (I) do not leave the trap. We shall see 
that this has important consequences for all other processes to be discussed. In the literature, the rate of process (I) 
is often denoted as L3. This rate has been calculated by Petrov [2]. The result is that L 3 is proportional to the a® 
and depends linearly on the average kinetic energy of the particles. 

Together with Process (I) is the reverse process where a dimer dissociates back to atoms through collisions with 
other atoms and with dimers, which we denote as 

A Ul) + D -> A Ui) + A t + A l Process (I'). (2) 

This will lead to an increase of the number of atom and a decrease of the number of shallow dimers (D) in the 
trap. Process (I') shares the same microscopic matrix element as (I). However, the density of states of the initial 
configuration of these two processes are entirely different. For the dissociation process (I'), it depends crucially on 
temperature and interaction parameters of the system. The population of shallow dimers in the trap is a result of 
the competing effects (I) and (I'). 

While both (I) and (I') change the numbers of atoms and dimers, none of these particles are lost from the trap 
during these processes, since £5 < VJj as mentioned before. Of course, shallow dimers can also dissociate through 
dimer-dimer collision, 

D + D ->• A t + A x + A t + A l Process (I"). (3) 

We shall ignore Process (I") since it has the same microscopic matrix element as its reverse process in which four 
atoms collide to form two shallow dimers, which is very small. While (I) and (I') do not cause particle loss from 
the trap, processes leading to formation of deep bound dimer state will. These processes arise from collision between 
shallow dimers and atoms, (referred to as Process (II)) or between shallow dimers with each other (referred to as 
Process (III))- Schematically, they are represented as 

A Hl) A Ul) + D * Process (II), (4) 

D + D^D*+D* Process (III). (5) 

where D* represents the deep bound state. Since the energy of the deep bound state is large and negative, the atoms 
and dimers in the final state in Process (II) and (III) will carry very large kinetic energy and will leave the trap. 
The rates of Processes (II) and (III), denoted as Li and L m , respectively, have also been calculated by Petrov et 
al. [13] . who found L2 cx aj 3 3 and L rn oc aj 2 5 . 

The Origin of the maximum in atom loss rate: The processes above show that the loss of atoms from a trap 
proceeds in two-steps. The first is to produce a density of shallow dimers, which is the combined effect of the three 
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body-recombination (I) and dimer disassociation (I'). Once the dimers are formed, Processes (II) and (III) will 
be activated and produce atoms and dimers with kinetic energies high enough for them to leave the trap. As one 
approaches the resonance from small k F a s , the rate of 3-body recombination (process (I)) increases as [H]. The 
population of the shallow dimer, and hence the atom loss rate, therefore increases with a s . However, as one gets 
closer to resonance, kpa s — > +00, the binding energy of the shallow dimer E(, = h 2 /ma 2 s decreases rapidly, making 
the dissociations process (I') more and more effective, caused by thermal effects at high temperatures, or by Pauli 
blocking effects in quantum degenerate regime |15j . (see later discussions). At some point, dimer dissociation (I') will 
overwhelm the 3-body recombination process (I), thereby quenches Processes (II) and (III) and reduces the loss rate. 



IV. THEORETICAL MODEL 



In this section, we express the processes discussed in Section III in terms of rate equations. Our basic assumption 
is that there is a well-defined degrees of freedom for the shallow dimers. We shall denote the number of atoms of up 
and down spins as n+ and ni, and the number of atoms and number of dimers are denoted as n a and n m respectively, 
and n a = + ri\,. The processes in Section III imply 



dn m (t) 
dt 



= [L 5 (n*(t)n\,(t) + n t (t)nl(t)) - qL 3 n a (t)n m (t)] - L 2 n a (t)n m (t) - 2L m n m (tf 



= [-13(71^(^(1) + n t (t)nl(t)) + qL 3 n a {t)n m (t)] - L 2 n t (t)n m (t) 
5r ^ ) = [-L 3 (^(t)n 4 (t) + n t (i)nJ(t)) + ? L3«a(*)rim(*)] - L 2ni (t)n m (t) 



(6) 
(7) 
(8) 



The first terms on the right hand side of Eqn.([6]) describes the production of dimers through three-body recombination 
(Process (I)) with rate L 3 . The second term describes the dissociation process (I') through atom-dimer collision 
(hence proportional to n a n m ). We have parametrized the rate of the process as L 3 q. As we shall see, q depends on 
temperature T, scattering length a s , and the total number of particles n = n a + 2n m ; 



q = q(T, a s ,n), 



(9) 



since the ratio between atom and shallow dimers in equilibrium depend on the total particle number n. (See 
discussions later in this section). That we group these two terms (L 3 and qL 3 ) together in a square bracket because 
they do not contribute to particle loss in the trap. (See discussions in the previous section). The third term in 
Eqn.Q describes the loss of shallow dimer in Process (II) through collisions between atom and shallow dimer that 
leads to production of deeply bound states. The rate of this process is L 2 . The last term in Eqn.(|6]) describes 
Process (III), which describes the loss of shallow dimers due to collisions between them to form deeply bound states. 
The rate of this process is L m . Eqn.([7]) and Q can be understood similarity, based on the picture that the dimers 
produced in Process (I) do not leave the system, (since the binding energy of the dimer Eb is much less than the 
typical trap depth shown in Table [l]). As a result, the shallow dimers produced continue to participate in the 3-body 
recombination process and dissociation process. 

With q being a function of n = n a + 2n m = iVf + n_i + 2n m , Eqn.(|6|, (|7j), and |8]) form a close set of equations for 
7tf, n.j,, and n m , and allow us to study the time evolution of these quantities. For the case of equal spin population, 
= ni = n a /2, Eqn.( 6|7p ) can be simplified to 



dn m (t) 

dt 
dn a (t) 
dt 



J 3 3 



n a (t) - qL 3 n a (t)n m (t) 



L2n a (t)n m (t) - 2L m n m (ty 



L 



3 3 



n 6 a (t) - 2qL 3 n a (t)n m (t) 



L 2 n a (t)n m (t). 



Note that from Eqn.( 10 ) and ( 11 ), the total number of particles n = n a + 2n m decreases as 

dn 



dt 



= -2L 2 n a (t)n m (t) - 2L m n m (t), 



(10) 
(11) 

(12) 



which is the statement that only Processes (II) and (III) lead to particle loss. 
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To make use these equations, we need to obtain expressions of the parameters L 3 ,L 2 , L m and the function q(T, a s ,n). 
The typical values of L3, L 2 and L m can be obtained from existing experiments. Our strategy is to fix their values 
at one scattering length a s and then use the scaling relation derived by Petrov et al. [T3] to get the value at other 
scattering length. For 6 Li, the values of L3 and L m are given in Jochim et al [TO]. At field B = 690G, we have 
L 3 (B = 690G) = 1 x 10- 25 cm 6 /s. Thus at a different magnetic field B, 



L 3 (B) 



a s (B) 



a s (B = 690G) 



L 3 (B = 690G) 



The value of L m is estimated to be 5 x 10 11 cm 3 /s in S. Jochim et al. [10] at field B 
calculation by Petrov et al. at other magnetic field B, the rate coefficient is given by 



(13) 



546G. According to the 



L m (B) 



a s {B) 



a s {B = 546G) 



-2.5 



L 2 (B = 546G). 



(14) 



We haven't been able to find a precise value for L 2 m existing experiments. As commented in ref.[S], the value of L 2 
cannot be safely estimated, but a reasonable value can be estimated to be L 2 (B = 690G) = 1 x 10 _13 cm 3 /s and thus 
at another magnetic field 



L 2 (B) 



a s (B) 



-, -3.3 



a s (B = 690G) 



L 2 (B = 690G). 



(15) 



Note that strictly speaking, the above formulae were derived for the regime kpa s < 1 |T3l HI]. We shall assume in 
our subsequent discussion that these expressions continue to give a reasonable approximation to actual rates in the 
region kpa s ~ 1. 

To determine q, we use the fact that in the absence of the L 2 and L m terms, the long time evolution of the equations 
above should establish chemical equilibrium between atoms and molecules. In that case, we find 



Aqn r , 



(16) 



where 



''a, eg 



limt-^oo n a (t) and similarly for n m ^ eq . Eqn.(16) shows that the quantity q is simply the ratio of atoms 
and molecules in an equilibrium mixture. To estimate n a ^ eq and n m ^ eq (and hence q), we use the simplest model of 
non-interacting mixture, with the approximate hamiltonian K = (h a — fi a n a ) + (h m — [i m n m)-, where h a and h m are 
the hamiltonians of the atoms and shallow dimers respectively, and \x a and [i m are their chemical potentials, (see also, 
for example, Kokkelmans et a/.[TT], as well as Chin and Grimmjl2j). The energy of the dimer will be denoted as Ef,. 
Equilibration between atoms and dimers are ensured through the relation of their chemical potentials 



2/i a — ~E b + fi m , 

and that the number of atoms n a and the number of dimers are constraint by the condition 



The explicit for of n a and 



d 3 k 



L rn.eq 



(27r) 3 Jk 2 /2m-n a ) 

d 3 k 1 

(27r) 3 (t 2 /4m- Pm ) _ 



(17) 
(18) 

(19) 
(20) 



- 1 

where =l/fcgT is the inverse temperature, and m is the mass of the atom. 

In vacuum, we have Ef, = h 2 /ma s . In a many-body system, as temperatures drops to quantum degenerate regime, 
pair formation is affected by the presence of a Fermi sea, which has taken up some momentum states needed for the 
pair wave function. To capture the effect of Pauli blocking on formation of bound state, we consider the analog of 
Cooper pairing in the presence of a Fermi sea at finite temperature. Denoting the quantum state of Cooper pair as 
|^) = J2k ^klk t ; k I), the Schrodinger equation of the pair is 



k' 



(21) 
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where "Jk is the amplitude for presence of a pair, 6k = ft 2 /2m, fl is the volume, /k is a Fermi function for the atom, 



1 



(22) 



and u is bare interaction parameter which is designed to reproduce the low energy scattering amplitude, and is 
related to the s-wave scattering length a s as 



to 
4Trh 2 a s 



k 

The reason for the power 2 in the Fermi exclusion factor 1 — /k' is because the scattered state consists of two particles 

■(Wk) 2 



k' and -k'. The solution of Eqn.plj) is 



TO 



-y 



1 



Eb — 2ek 2ek 



(24) 



Eqn.(24| gives Eb as a function of T, a s and [i a . It is easy to see that in in the non-degenerate limit, e M °/ T — » 0, we 
have /k — > 0, and Eqn.(24) reduces to the equation for bound state energy in vacuum, which gives the usual result 
E b = h 2 /(ma 2 s ). 



Eq.(18) and (|17|) now imply 



d 3 k 



d 3 k 



(27r) 3 e (i 2 /2m- (Ia ) J (2ir) 3 e (k 2 /4m-2tJ, a -E b ) _ ^ ' 



(25) 



Since Ef, is a function of T, a s and /U a , Ej, = Eb(a s ,T, (i a ). Eq.(25l gives /i a as a function of n, T, and a s . Once 
Ma = fJ-a{n,T, a s ) is determined, we can then calculate n a and n m from Eqn.(19) and (20), and obtain q Eqn(16), 



g(T, a s , n) = n 2 a / (4n m , eg ). 



(26) 



Before proceeding, let us comment on the assumptions of Eqns.( 10|11 1. (i) First of all, we have assumed a well 
defined degrees of freedom for the shallow dimers, which is only valid for small kpa s or high temperature. As a 
result, our equations will not be accurate very close to the resonance where the dimer degrees of freedom becomes 
less well defined due to many-body effects, (ii) The most important effect that is not captured by Eqns.( T0|11 1 is 
heating. In principle, there should be an equation of the form dT/dt — F(T, n a ,n m ; a s ), where F describes the effect 
of re-thcrmalization of the energy released form various processes. We have not attempted to construct this equation. 
The viewpoint we take is that since the physics of the two-step process has already provided an explanation for the 
loss maximum qualitatively, (see Section [V]), it is useful to first find out how well these equations account for current 
experiments quantitatively so as to determine the validity and usefulness of these rate equations. The more elaborate 
effects of heating will be explored elsewhere, (iii) Throughout our discussion, we shall replace the quantum gas in a 
harmonic trap with a non-uniform density profile by one in a square box with uniform density. The replacement is 
mainly for simplicity. While one can perform a full calculation using local density approximation (LDA), for the level 
of accuracy of our comparison, we believe our simple replacement is sufficient. 



V. COMPARISON WITH EXPERIMENTS 



Essentially, all experiments on atom loss are performed in the following way. One starts with a sample with fixed 
number of atoms n and no dimers at an initial magnetic field Bi„i corresponding to a small and positive scattering 
length a" 11 . One then jump to the magnetic field Bfi n , corresponding to the scattering length a s of interest and wait 
for a time interval At. Within this time interval, atoms are converted to shallow dimers as well as escaping from the 
trap due to formation of deep bound states. At the end of this interval, the system has n a atoms and n m shallow 
dimers. The number of atoms n a is the quantity of interest. To image the number of remaining atoms after the 
interval At, one pulls back the magnetic field B from Bji n to B rneasure , typically corresponds to a scattering length 
that is small and positive a ™ easure . The rate of the pull back is such that all the shallow dimers formed at magnetic 
field Bfi n are converted to tightly bound state of size a ™ easure and therefore will not be counted by the imaging 
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FIG. 1: The calculated atom remaining in the trap for the MIT 2002 experiment. The initial density is given by n a (t = 0) = 
3 x 10 13 cm -3 and n m (t = 0) = 0. The corresponding function q(T,a B ) for the appropriate density is shown as well. 



process for atoms. The number of atoms counted at B measure (or at a ™ easure ) therefore gives the number of atoms 
n a at the end of the time interval At. 

We have applied the rate equations in Sec IV to different experiments listed in Table [I] The results and the 
parameters used in our calculations are summarized below: 



A. MIT/2002 

In this experiment [6j, the initial ensemble consists of 3 x 10 5 6 Li atoms with peak density 3 x 10 13 cm~ 3 . The 
temperature where the final measurement is made is around T = 22/iK, slightly higher than the Fermi temperature 
Tp = 21fiK. The trapping parameters are r = 12Hz, Z = 200Hz and the trap depth Vq = Ylh^iK. The magnetic 
field is turned on within 4ms to -B/m and one waits for 50 ms or 500ms at B = Bfi n , then the magnetic field is 
switched off within 100/Lts and the cloud is probed by absorption imaging. The experimental findings are that for 
incoherent two-component Fermi gas, there is a strong loss around magnetic field B = 680G. Also, for B > 680G 
and close to the unitarity Bq = 834G, the loss decreases and saturates. 

To describe the experiment, we choose as our initial conditions for the rate equations n = n a (t — 0) = 3 x 10 13 cm -3 
and n m (t = 0) = 0, corresponding to the experiment. We run our rate equations for Ai = 0.045s. We have chosen 
the following parameters for the experiment. L 3 (B = 690G) = 1 x 10~ 24 cm 6 /s, L 2 (B — 690G) = 1 x 10 _13 cm 3 /s 
and L m (B = 546G) = 5 x 10~ 11 cm 3 /s. The calculated fraction of atoms remained is shown in Figure [l] We note 
that the atom loss behavior in this experiment is quite insensitive to the value of L rn . 



B. Innsbruck/2003 

This set of experiments 10J are described in detail in S. Jochim's thesis[TD]. One usually starts with about two 
million atoms in the lowest two hypcrfine-Zeeman states of the 6 Li atoms. The samples are cooled to three different 
temperatures: 22/iK, 30/iK and 60/iK, at a magnetic field 300G where the scattering length is large and negative. As 
a result, there are no shallow dimers in the initial state. One then ramps the system to close to Feshbach Resonance 
within 50ms. After waiting for 5 ~ 7s, the magnetic field is ramped back to zero, at which point the number of 
atoms are measured. The atom loss is found to have a maximum at 636G. To describe the experiment, we run 
the equation for At = 5s, with n — n a (t = 0) = 3.8 x 10 13 cto -3 , n m (t = 0) = estimated from ref. [10] . and 
L 3 {B = 690G) = 8 x 10~ 25 cm 6 /s, L 2 (B = 690G) = 1 x 10- 14 cm 3 /s and L m {B = 546G) = 5 x 10- u cm 3 /s'. The 
calculated number of atoms remained is shown in Figure [2] 




FIG. 2: The calculated atom remaining in the trap for the Innsbruck 2003 experiment. The initial density is given by 
n a (t = 0) = 3.8 x 10 13 cm -3 and n m (t = 0) = 0. The corresponding function q(T,a s ) for the appropriate density is shown as 
well. 



C. JILA/2004 

The starting point of this experiment [7] is an ensemble of 40 K atoms in the hyperfine-Zeeman states 1 9/2, —9/2) 
and |9/2, -7/2), with temperature T = 70nK and T/T F = 0.22. The peak density of the system is 1.5 x 10 13 cm" 3 . 
In this experiment, one uses the radio-frequency spectroscopic to disassociate the molecules in the states |9/2, —9/2) 
and |9/2,— 7/2) to atoms in the hyperfine-Zeeman states 1 9/2, — 9/2) and 1 9/2 , — 5/2). By measuring the atoms in 
the state 1 9/2, —5/2), one infers the molecule numbers of the system. 

The atom loss shows the expected non-monotonic behavior for both spin component. To study this experimental 
situation, we choose the temperature of the system to be 0.67fiK, and run our rate equation for At = 95ms, which 
is the holding time in the experiment. The maximum loss occurs at around 201. 5G. The following parameters 
are assumed for the experiment. L 3 (B = 200G) = 1.5 x 10~ 25 cm 6 /s, L 2 (B = 200G) = 2 x 10~ 12 cm 3 /s and 
L m (B = 200G) = 1 x 10~ n c7n 3 /s. We have taken n = n a (t = 0) = 1.5 x 10 13 cto~ 3 , n m (t = 0) = 0. The calculated 
atom loss rate is shown in Figure [3j 




FIG. 3: The calculated decay rate for the JILA 2004 experiment. The initial density is given by n a (t = 0) = 1.5 x 10 cm 
and rim.it — 0) = 0. The corresponding function q(T, a s ) for the appropriate density is shown as well. 



D. MIT/2009 

In this experiment [5j, the Fermi gas is prepared initially at a field B around 600G. The central density of the 
trapped gas is n(r = 0) = 0.69 x 10 13 cm~ 3 . The trap frequencies are v x = v v = 300Hz and v z = 70Hz . In this 
experiment, the field is ramped down within 4.6ms to the desired final field. It is noted that for kpa s > 1.8, there 
is approximately 25% percent molecule population after the ramp, independent of temperature (within the range 
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FIG. 4: Calculated decay rate for the MIT 2009 experiment. The initial atomic density is about n a (t = 0) = 0.69 x 10 13 cm -3 . 
The three curves correspond to different choice of the positions in the trap. Blue: r* = ^/2/5Rtf', purple: r* = ^1/2Rtf 
and brown: r* = ^/3/5-Rtf- The corresponding curve for q(T, o s ) is shown as well. 



considered). A maximum in atom loss similar to the observations of other groups is found, but occurs around 780G. 
From the information in Ref. [5], we take T = 0.3(iK. To study this experiment using our rate equation, we choose 
L 3 (B = 690G) = 1 x 10- 29 cm 6 /s, L 2 (B = 690G) = 1 x 10- 13 cm 3 /s and L m (B = 546G) = 1 x 10- n cm 3 /s. We run 
our equation for At = 2ms, since it was remarked in Ref. [5] that the rate are "measured within the first 2 ms". We 
have run our equations with three different initial densities n — n a (t — 0) chosen as follows. Since we have replace 
the trapped gas in a harmonic well by one in a cubic box, we take n as some average density of the density in the 
harmonic trap, i.e. setting n — titf{t*), where titf{f) is the Thomas- Fermi density profile of the Fermi gas in the 
harmonic trap, and r* is some radius less than the Thomas-Fermi radius Rtf ( r * < Rtf)- The calculated atom loss 
rate is shown in Figure [4] The three curves in Figure [4] are different choices of n corresponding to different choses of 
r* . The color scheme in Figure [4] is that blue, purple, and brown correspond to r* — ^/2/5Rtf, r* — \J\I2Rtf, and 
r* = y^3/5R TF , respectively. 



E. ENS/2004 

This experiment [5] starts with a gas of 7 x 10 4 6 Li atoms in anisotropic trap with ui x — 2tt x 0.78kHz, 
ojy = 2tt x 2.1kHz and ui z — 2ir x 2.25kHz. The temperature of the system is estimated to be T = 2.4/iK, while 
the Fermi temperature Tp sa 6/iK. The system is prepared at low magnetic field and then evaporative cooling is 
performed at field B — 320G, where the scattering length a s = — 8nm. In 10ms, the magnetic field is ramped to 
anywhere between 600 and 850G and time-of-flight expansions are taken and number of atoms is counted. It is 
found that the loss rate has a maximal around B = 720G. Interestingly, together with this loss maximum, the sign 
of interaction energy changes at the same point. 

The following numbers are used to fit the experiment. The central density is given by n(r = 0) = 3.5 x 10 13 cm -3 . 
L 3 (B = 690G) = 1.5 x 10~ 26 cm 6 /s, L 2 (B = 690G) = 4 x 10- 13 cm 3 /s and L m (B = 546G) = 1 x 10- 12 cm 3 /s. We 
run the equation for At = 1.5s. We have not been able to find the holding time in ref. [8]. However, it is remarked 
in the cited reference that the life time of the gases ranges from 100ms to few seconds. Our choice of 1.5s is a rough 
estimated of the holding time. We have checked that a different choice of reasonable holding time does not change 
the loss behavior, provided that we modify appropriately other parameters in the calculation. The calculated atom 
loss rate is shown in Figure [5j 

Finally, we would like to point out that in all our studies, the emergence of a maximum atom loss as a function of 
a s is a robust phenomena. Modest changes of input parameters, as well as the rates L 2l L 3 and L m do not change 
our results. 
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FIG. 5: The calculated atom remaining in the trap for the ENS 2004 experiment. The initial density is given by n a (t = 0) = 
3.5 x 10 13 cm~ 3 and n m (t = 0) = 0. The corresponding function q(T,a s ) for appropriate density is given by as well. 



VI. CONCLUSIONS 

In this work, we investigate the origin of the maximum of atom loss as a function of scattering length observed in 
many experiments that are performed over a wide range of temperature, trap depths, and particle numbers. We find 
that this is a result of a two-step process. The first is the production of a population of shallow dimers which remains 
in the trap. The second is the conversion of these shallow dimers into deep bound states through collision processes. 
The second step, which causes particle loss from the trap, depends on the number of shallow dimers generated in the 
first step. The maximum of atom loss is caused by the variation of the number of shallow dimers as a function of a s . 
This number is the result of the competition between 3-body recombination and the dimer dissociation processes. 
While the former increases the dimer population at a rate rises as af, the later reduces it and becomes more and 
more effective as one approaches the resonance, since the binding energy of the dimer decreases rapidly. We have 
cast these processes in a set of simple rate equation, and show that they can account for the atom loss observed in 
all current experiments. 

Apart from these agreements, we note from Table [i] that E^/V ~ 0.1 for all current experiments. We believe 
this is not an accident. Typically, because of the evaporation process, the temperature of the system is lowered if 
the trapped depth is lowered. Let a* be the scattering length of the maximum atom loss at temperature T. By 
raising the temperature (caused by a higher depth V ), one makes the dissociation process more effective. As a 
result, one will have to go to a small scattering length (hence larger binding energy E£ = h 2 /ma 2 ) to achieve the 
same ratio between the rate three body recombination and the rate of dissociation. This shows that the observed 
maximum of atom loss is not purely a function of kpa Sl but depends on external parameters like trap depth 
and temperature. In other words, even though experiments are performed with the same atom density kp, the 
location of the maximum loss will be different for different runs with different trap depth V , and the result would 
appear irreproducible if variations of extrinsic factors are not fully taken into account. In any case, our results show 
that the maximum of the atom loss can not be used as a tool to determine the nature of the ground state of the system. 

This work is supported by NSF Grant DMR-0907366 and by DARPA under the Army Research Office Grant Nos. 
W911NF-07-1-0464, W911NF0710576. 
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